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ABSTRACT 

We present multi-band optical and near-infrared observations of 19 short 7-ray burst (GRB) host galaxies, 
aimed at measuring their stellar masses and population ages. The goals of this study are to evaluate whether 
short GRBs track the stellar mass distribution of galaxies, to investigate the progenitor delay time distribution, 
and to explore any connection between long and short GRB progenitors. Using single stellar population models 
we infer masses of log(M„/M Q ) ps 8.8- 11.6, with a median of (log(M„/M Q )) ps 10.1, and population ages 
of t* fa 0.03-4.4 Gyr with a median of (r*) ps 0.3 Gyr. We further infer maximal masses of log(M*/M & ) ps 
9.7- 1 1 .9 by assuming stellar population ages equal to the age of the universe at each host's redshift. Comparing 
the distribution of stellar masses to the general galaxy mass function we find that short GRBs track the cosmic 
stellar mass distribution only if the late-type hosts generally have maximal masses. However, there is an 
apparent dearth of early-type hosts compared to the equal contribution of early- and late-type galaxies to the 
cosmic stellar mass budget. These results suggest that stellar mass may not be the sole parameter controlling 
the short GRB rate, and raise the possibility of a two-component model with both mass and star formation 
playing a role (reminiscent of the case for Type la supernovae). If short GRBs in late-type galaxies indeed 
track the star formation activity, the resulting typical delay time is ~ 0.2 Gyr, while those in early-type hosts 
have a typical delay of ~ 3 Gyr. Using the same stellar population models we fit the broad-band photometry 
for 22 long GRB host galaxies in a similar redshift range and find that they have significantly lower masses 
and younger population ages, with (log(M»/M©)) « 9.1 and (r*) ps 0.06 Gyr, respectively; their maximal 
masses are similarly lower, {log(M*/M Q )} ps 9.6, and as expected do not track the galaxy mass function. 
Most importantly, the two GRB host populations remain distinct even if we consider only the star-forming 
hosts of short GRBs, supporting our previous findings (based on star formation rates and metallicities) that the 
progenitors of long GRBs and short GRBs in late-type galaxies are distinct. Given the much younger stellar 
populations of long GRB hosts (and hence of long GRB progenitors), and the substantial differences in host 
properties, we caution against the use of Type I and II designations for GRBs since this may erroneously imply 
that all GRBs which track star formation activity share the same massive star progenitors. 

Subject headings: gamma-rays:bursts 
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1. INTRODUCTION 

The progenitors of the short-duration 7-ray bursts (GRBs) 
remain unidentified at the present. However, in recent years 
important constraints have been placed on their nature based 
on observations of the prompt emission, afterglows, and host 
galaxies. In p articular, the discovery of short GRBs in ellip- 
tical galaxies dBerger et al.ll2005t iBloom et al .1120061) demon- 
strated that at least some of the progenitors belong to an old 
stellar population (> 1 Gyr). On the other hand, the bulk of 
the identified host p opulation exhibits star formation activ- 
ity dBerger et al.ll2007t) . albeit with lower star formatio n rates 
and h igher metallicities than the hosts of long GRBs (Berger 
2009). On a sub-galactic scale, the locations of shor t GRBs do 
not coincide with regions of active star formati on dFox et al.1 
120051; iFong et al.ll2010t iRowlinson et~ai1l2010bl) . and instead 
trace the rest-frame optical light distribution, indicative of a 
connec tion with older stellar population s compared to long 
GRBs (lFongetalJl20T0t lBergerll20T0bl) . Similarly, the off- 
sets of short GRBs are s ubstantially larger than those of long 
GRBs (Fong et al. 2010), with recent potential evidence for a 
population of shor t bursts with of fsets of tens of kpc, possibly 
due to natal kicks dBergerCOlOal) . 

The afterglow and prompt emission properties likewise 
point to a distinct origin for the two GRB types. Short 

1 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138 



GRBs have significantly fa inter op tical a nd X-ray afterglows 
than those of long GRBs (iNakarl 120071; iGehrels et al l 2008 
i Kann et al.ll2008HBerger et alj|2009t INvsewander et alj|2009 



lBergerll2010al) . altho ugh they exhibit similar flux ratios be - 
tween the two bands (Nysewa nder et al .ll2009b iBer gerll20 1 Oah . 
Their isotropic-equivalent 7-ray energies, £ 7i i s0 ~ 10 49 — 10 52 
erg, a r e also generally lowe r than for long GR Bs ( Berg eret al.| 
1 20071: IGehrels etalJ 120081: INvsewander eTail l200ft iBerger 
2010a). These results point to an overall lowe r energy budget, 
and potentially low er circumburst densities dSoderberg et ID 
l2006t feerger 2010a). Taken in conjunction, the emission and 
environmental properties point to a distinct origin for short 
and long GRBs, and indicate that the short GRB progenitors 
are not dominated by a young population (e.g., massive stars 
or young magnetars). 

Despite this crucial insight three basic and inter-related 
questions remain open at the present, whose answers can shed 
light on the nature of the progenitors: 

• Does the short GRB rate depend solely on the stellar 
masses of the host galaxies, or does star formation play 
a role? 

• What is the delay time distribution of short GRBs? 

• Is there an overlap between long and short GRBs, such 
that duration is not uniquely mapped to progenitor? 
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In the context of the first two questions, the popul ar compact 
objec t binary progenitors (NS-NS or NS-BH; lEichler et all 
Il989t iPaczvnskil [19911: iNaravan et ail 1 19921) are expected to 
track stellar mass if the merger timescales are long (> Gyr). 
On the oth er hand, if a substanti al rapid merger channel ex- 
ists (e.g., iBelczynski et al.N2006l) . then star formation activ- 
ity will also impact the short GRB rate. Similarly, the dis- 
tribution of merger t imescales will impact the mix of early- 
and late-type hosts (Zheng & Ramirez-Ruiz 2007). In par- 
ticular, if mass is the dominant factor we expect a roughly 
equal number of short GRBs in early- and late-type galax- 
ies 2 since at the relevant redshift range each galaxy population 
type acc ounts for about h alf of the cosmic stellar mass bud- 
get (e.g.. lllbert et al.ll2010b ; a connection with star formation 
will serve to increase the late-type fraction. The same trends 
will be apparent for any progenitor system with a broad range 
of delay time s, such as accre tion-induced collapse (AIC) of 
neutron stars dOin et al.l lT998). and delayed magnetar forma- 
tion through binary white dwarf merge rs or white dwarf AIC 
(iLevan et al.ll2006t iMetzger et al.ll2008l) . 

Related to the questions of delay times and a connection to 
stellar mass or star formation is whether durations and pro- 
genitors are directly linked, i.e., whether some short GRBs 
arise from massive star progenitors l ike long GRBs (e.g., 
I VirgiU et al.l l2009t lLazzati et all 1201 Ot) . If such a subset of 
events exists, it will be manifested in a population of short 
GRBs with negligible delay times of ~ 10 Myr, or equiv- 
alently hosts with very young stellar populations and vig- 
orous star formation activity. Suc h stellar populations are 
clearly present in lo ng GRB hosts (IChristensen et al.l)2004j: 
iLevesque et al.ll2010al). In this context. Khang et all (2007) 
and |Zhang et al.l d2009l) proposed a re-classification of long 
and short GRBs into Type I and II events associated with old 
and you ng stellar pop ulations, respectively. Building on this 
concept, Virgili et al.l (120091) recently carried out a joint anal- 
ysis of the short GRB luminosity -redshift distribution and the 
BATSE logN-logS distribution and concluded that ~ 90% 
of the bursts classified as short have negligible delay times 
and therefore directly track the cosmic star formation history. 
These authors therefore conclude that most short GRBs have 
massive star progenitors similar to those of long GRBs. 

Here we take a different approach to exploring the relation 
between short GRBs and stellar mass, the short GRB delay 
time distribution, the possibility that short GRBs also track 
star formation, and whether such a connection indicates a re- 
lation to long GRB progenitors. Our study relies on multi- 
band optical and near-infrared observations of short GRB host 
galaxies from which we extract the relevant stellar population 
properties: mass and age. We confront the resulting mass 
distribution with the general galaxy stellar mass function to 
assess the relation to mass for the population as a whole, 
and separately for the early- and late-type hosts. We also 
use the inferred stellar population ages to place constraints 
on the typical delay times in early- and late-type hosts. Fur- 
thermore, we uniformly re-analyze the optical/near-IR spec- 
tral energy distributions of long GRB hosts to allow a direct 
comparison with the inferred properties of short GRB hosts. 
This allows for a comparison of the stellar populations in 
the two host samples with greater diagnostic power than the 
use of secondary indicators such as the luminosity-redshift or 

- This assumes that there are no systematic differences between short 
GRBs in early- and late-type galaxies that would influence their detectability, 
such as the typical ambient density or energy release. 



logiV-logS distributions. 

The plan of the paper is as follows. The sample properties 
and optical/near-IR observations are summarized in $2] The 
stellar population models and fitting procedure are described 
in $3] In $4] we describe the resulting distributions of stellar 
masses and ages for a range of models and assumptions. We 
confront the mass distribution of short GRB hosts with the 
galaxy mass function in jj5.ll and compare the stellar popula- 
tion properties of long and short GRB hosts in 35.21 We draw 
conclusions about the delay times of short GRB progenitors 
in the context their relation to mass and star formation in 35.31 
and summarize our key findings in ^6] 

2. SHORT GRB SAMPLE AND OBSERVATIONS 

The sample of short GRB hosts described in this pa- 
per includes 19 galaxies, of which 9 have been iden- 
tified in coincidence with optical afterglow positions 



tinea in co i nciden c e witn optical aiterglow po 
dBerger et alJ 120051: iFox et alj 120051: iHiorth et alj 
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Covino et al . 2006; Soder berg et al. 2006; Be rger et alJ 
D'Avanzo etail 120091: iGraham et all l2009t iMcBreen et alJ 
l2010tlRowlinson et a l. 2010a), 8 have been identified in coin- 
cidence with Swift / X-my Telescope (X RT) error circles (e.g., 
iBloom et alJl2006l: iBerger et alJl2007h . and 2 (GRBs 070809 
and 090515) have been identified as putative hosts based on 
chance coincidence probabilities for bur sts with optica l after- 
glows and no coincident host galaxies (Berger 2010a). Two 
of the hosts (051210 and 070729) do not have spectroscopic 
redshift measurements, and we infer photometric redshifts as 
part of the stellar population analysis presented here (§0. In 
the analysis below we address the properties and implications 
of the full sample, as well as the subset of secure host associ- 
ations (i.e., those found in coincidence with optical afterglow 
positions). The properties of the bursts are summarized in Ta- 
ble Q] We note that a few additional secure hosts have photo- 
metric measurements in only 1-2 optical filters and no spec- 
troscopic redshifts (e.g., GRBs 060121 and 060313); these do 
not allow us to infer the required photometric redshifts, stellar 
masses, and ages with any confidence. 

2.1. Optical and Near-IR Observations 

We obtained multi-band optical and near-IR observations 
that cover rest wavelengths across the Balmer break and 
in the near-IR. This allows us to robustly assess the stel- 
lar population ages and the stellar masses. Optical observa- 
tions were obtained with the Low Dispersion Survey Spec- 
trograph (LDSS3) on the Magellan/Clay 6.5-m tel escope, and 
the G emini Multi-Object Spectrograph (GMOS; iHook et al] 
2004) on the Gemini North and South 8-m telescopes. Near- 
IR observations were performed with the Persson's Auxiliary 
Nasmyth Infrared Camera (PANIC) on the Magellan/Baade 
6.5-m telescope. 

Reduction and photometry of the Magellan observations 
were performed using standard routines in IRAF, including 
dark frame subtraction and fringe correction for the PANIC 
data. The Gemini data were processed using the gemini 
package in IRAF. The resulting photometric measurements 
(in AB magnitudes and corrected for Galactic extinction us- 
ing Schlegel et al. 1998) are listed in Table Q] and the spec- 
tral energy distributions (SEDs) are shown in Figure [T] As 
noted in Table Q] our host observations ar e supplemented by 
data from the litera ture for G RBs 050709 dHiorth et 111120051: 
ICovino et al.ll20 06). 050724 dGorosabel et al.ll2006l). 05 1221 
dSoderberget all 120061). 061006 dD'Avanzo et al.ll2009l). and 
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070714b (Grah am etal]l2009l) : for the host of GRB 050509b 
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we used photometric measurements from the SDSS and 
2MASS catalogs. 

To convert fluxes to luminosities we use the standard cos- 
mological parameters: Ho = 70 km s" 1 Mpc" 1 , Sl m = 0.27, and 
n A = 0.73. 

3. STELLAR POPULATION MODELING 

To determine the stellar masses and population ages of the 
short GRB host galaxies we fit the optical/near-I R SEDs with 
the stellar population models of Maraston (2005). These mod- 
els provide the spectral luminosity as a function of wavelength 
for an equivalent of 1 M Q . We use the single age simple stel- 
lar population (SSP) red horizontal branch morphology mod- 
els with a Salpeter initial mass function. We further use So- 
lar metallici ty as generally appropriate for short GRB hosts 
(Berger 2009); we note that changing the metallicity to 0.5 
Zq results in a systematic decrease in the inferred masses and 
ages by at most ~ 0.1 dex. 

To fit the observed SEDs with the SSP models, which exist 
for a discrete set of ages, we create an interpolated grid with 
several hundred logarithmically-spaced ages between 1 Myr 
and the age of the universe at each host's redshift. We then 
transform the model spectral luminosities and wavelengths to 
the observer frame and convolve with the appropriate filter 
response functions to create model magnitudes. The best-fit 
model for each host is determined using \ 2 minimization with 
the age of the population as the free parameter and the best-fit 
mass for each age (i.e., the normalization) determined by the 
equation: 



Mm = 



Xw=l ^V.model.i X ^V,data,i/°' w dataj 



V w F 2 

Z— d=l v. model, i 



(i) 



i/, data. i 



where F„ jm0< j B i i i are the model fluxes, F Wj data,i and cr„,data,i are 
the observed fluxes and uncertainties, respectively (Table [TJ, 
and N is the number of data points. The resulting best-fit sin- 
gle age SSP models are shown in Figure [1] and the inferred 
masses and ages are listed in Table [2] 

The SSP models provide a robust estimate of the stellar 
masses and population ages for the early-type host galaxies, 
since they are dominated by old stellar populations, which 
exhibit only a mild variation in the mass-to-light ratio as a 
function of age. However, for the late-type hosts, which have 
on-going star formation, the SSP fits essentially provide light- 
weighted (rather than mass-weighted) values, and they can 
therefore severely under-estimate the total stellar mass and 
the mass-weighted age of the stellar population. These are the 
quantities we are interested in for a determination of whether 
short GRBs track stellar mass. We therefore use two ap- 
proaches to assess the impact of older stellar populations for 
these hosts. 

First, we determine the maximum possible masses by as- 
suming that the mass of each host galaxy is dominated by a 
stellar population with the age of the universe at its redshift. 
Using the inferr ed (maximal) mass-to-light ratios from the 
Maraston (2005) models we use the available /f-band fluxes 
to extract the maximal masses. The resulting values are listed 
in Table [2] and provide an upper bound on the mass of each 
galaxy. An example of this approach compared to the single 
age SSP model is shown in Figure|2] 

Second, to provide a more realistic estimate we use a com- 
bined young+old SSP model to fit the entire SED of each late- 
type host galaxy and hence to determine the breakdown of 
masses in the young and old populations. The old population 



is fixed at the age of the universe at each host redshift, and 
we fit for the age of the young population and the mass ra- 
tio of the two populations; the old population mass is given 
by Equation Q] and the young population mass then trivially 
follows from the best-fit ratio. The results are summarized 
in Table |2] and an example of this approach compared to the 
single age SSP model is shown in Figure [2] We note that for 
some of the hosts a wide range of solutions is possible, in- 
cluding negligible mass in the old population (i.e., a solution 
identical to the original single age SSP fit). For these hosts we 
list the maximum possible mass in the old component within 
la- of the best fit. As expected, the young+old models lead 
to total masses intermediate between the single age SSP and 
the maximal masses. The addition of the old population also 
leads to systematically younger ages for the young population 
relative to the best-fit single age SSP values. 

With the exception of the host of GRB 071227 we find that 
host extinction is not required for a satisfactory fit of the ob- 
served SEDs. The addition of host extinction will generally 
result in systematically younger ages and somewhat higher 
masses for the young stellar population due to the increase in 
UV luminosity. On the other hand, extinction does not signif- 
icantly affect the total stellar mass since it is determined by 
the near-IR luminosity. We find that the addition of Ay w 0.5 
mag in the rest-frame of the host tends to change the SSP 
ages and masses (or the young population ages and masses in 
the young+old models) by at most -0.3 dex and +0.1 dex, re- 
spectively. Larger extinction values significantly degrade the 
goodness of fit. 

Finally, to uniformly compare the stellar masses and ages of 
short and l ong GRB hosts , we fit the SEDs of the long GRB 
hosts from Savaglio et al.] d2009l) in a similar redshift range 
(z ~ 0.03- 1.6) with the same stellar population models and 
techniques used here, but with a sub-solar metallicity. We 
derive single age SSP values, and maximal masses using the 
A'-band photometry alone. The results are listed in Table [3] 
As in the case of the short GRB hosts, we do not include the 
effects of extinction, unless strictly required by the goodness 
offit (see Table[3]). 

4. HOST GALAXY STELLAR MASSES AND AGES 

The distribution of stellar masses derived from the single 
age SSP fits is shown in Figure [3] We find that the masses 
span three orders of magnitude, Mssp ~ 6 x 10 s -4 x 10 1 1 M Q , 
with a median value of (Mssp) w 1.3 X 10 10 M Q . Dividing the 
sample into early- and late-type host galaxies we find that the 
former span the range M SSP sa (2-40) x 10 10 M , while the 
latter have much lower masses of Mssp ~ (0.06-2) x 10 10 
Mq. The clear distinction between the two samples partially 
reflects the bias of single age SSP models, which for the late- 
type hosts are dominated by the young stellar population and 
hence under-estimate the contribution of any older stellar pop- 
ulations. If we restrict the sample to the host galaxies identi- 
fied in coincidence with optical afterglows we find a median 
mass of (Mssp) ~ 5 x 10 9 M Q . This value is lower than for 
the full sample because 3 of the 5 early-type hosts, which 
have the largest SSP masses, are associated with XRT po- 
sitions (GRB 050509b) or inferred from chance coincidence 
arguments (GRBs 070809 and 090515). 

The differences between the masses of the early- and late- 
type hosts are less severe in the case of the maximal masses 
(Figure @J. For the population as a whole we find a reduced 
range of M Ma x ~6x 10 9 - 8 x 10 11 M Q . The median mass is 
(Mmslx) ~ 1 x 10 11 M Q , about an order of magnitude larger 
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than for the single age SSP masses, and only slightly larger 
than the stellar mass of the Milky Way. As expected, the ratios 
of maximal to SSP masses for the early-type hosts are modest, 
Mmux/Mssp ~ 2-8, since these hosts are already dominated 
by old stellar populations. However, for the late-type hosts the 
corrections are significant, Mmsx/A^ssp ~ 5-60, with a me- 
dian ratio of about an order of magnitude. If we again restrict 
the sample to host galaxies identified in coincidence with op- 
tical afterglows we find a median mass of (M Max ) w 3 x 10 10 
M©, lower than for the full sample due to the preferential re- 
jection of early-type hosts. 

While the maximal masses provide a robust upper bound, 
they are likely to over-estimate the true stellar mass, partic- 
ularly for the early-type hosts which do not suffer from the 
young-population bias present in the late-type hosts. If we in- 
stead use the young+old models for the late-type hosts, and 
the single age SSP values for the early-type hosts, we find 
masses of M w 2 x 10 9 -4 x 10" M , with a median of 
(M) « 5 x 10 10 M Q . The cumulative distributions of stellar 
masses for the host population as a whole, and for the subsets 
of early- and late-type galaxies are shown in Figure [5] The 
Figure clearly highlights the larger masses of the early-type 
hosts; even when comparing SSP masses for the early-type 
hosts with maximal masses for the late-type hosts the median 
for the former is larger by about a factor of 5. 

The range of inferred stellar masses does not appear to de- 
pend on redshift. As shown in Figure [6] the stellar masses are 
nearly uniform with redshift, both for the single age SSP and 
maximal values. We do find a preference for lower redshift 
among the subset of early-type hosts, and since these galax- 
ies are more massive, they lead to somewhat higher masses at 
lower redshift (z ~ 0.4 vs. z ~ 0.9). A similar conclusion is 
reached even if we restrict the sample to only the hosts iden- 
tified in coincidence with optical afterglows. 

We next turn to the distribution of stellar population ages. 
These values are only available for the single age SSP mod- 
els since for the maximal and young+old models we assume 
a population with the age of the universe at each host redshift. 
The distribution of ages is shown in Figure [7] with the values 
ranging from about 30 Myr to 4.4 Gyr. The median age is 
(tssp) ~ 0.3 Gyr for the full sample, with (tssp) ~ 0.25 Gyr 
for the subset of late-type hosts and (tssp) ~ 3 Gyr for the 
subset of early-type hosts. The median value for the subset of 
hosts identified in coincidence with optical afterglows is iden- 
tical to the median for the late-type hosts. As indicated above, 
the young+old models result in somewhat younger ages for 
the young stellar population since some of the near-IR flux 
is accounted for by the forced old population. The median 
age for the young population in this model is about 0. 1 Gyr 
(TableEJ. 

Finally, in Figure [8] we plot age as a function of redshift 
for the single age SSP models. For the full sample we find 
a significant negative correlation between the two quantities, 
with a Kendall's r value of -0.57, corresponding to a null 
hypothesis (no-correlation) probability of only m 0.01. How- 
ever, the correlation appears to be due mainly to hosts identi- 
fied in coincidence with XRT error circles and through chance 
coincidence probabilities. If we restrict the sample to hosts 
identified in coincidence with optical afterglows, then no sig- 
nificant correlation is found (Kendall's r value of -0.33 with a 
no-correlation probability of 0.26). We note that for the max- 
imal and young+old models we cannot associate an age with 
the old stellar population due to our conservative assumption 
that it is equivalent to the age of the universe. 



4.1. Long GRB Host Galaxies 

For the long GRB host comparison sample we derive single 
age SSP parameters, as well as maximal masses. For the SSP 
model we find a range of M S sp «6x 10 6 -2 x 10 10 Mm, with 
a median value of (M S sp) w 1.3 x 10 9 M Q (Figure |3). The 
maximal masses span M Ma x ~ 9 x 10 7 -9 x 10 10 M Q , with a 
median value of (M Ma x) ~ 4.0 x 10 9 M Q (FigurelU). Our me- 
dian SSP and maximal values bracket the median masses of 
2 x 10 9 M W and 1.8 x 10 9 M fl found bv lSavaglio etafl d2009b 
and iLevesque et alJ (1201 Obh . respectively. We also note that 
the median ratio of MMax/^ssp ~ 3 for the long GRB hosts 
is more modest than for the late-type short GRB hosts (see 
above), indicating that long GRB hosts do not generally har- 
bor massive old stellar populations. 

As in the case of the short GRB hosts, we find no clear cor- 
relation between redshift and stellar mass, although the low- 
est redshift hosts (z < 0.25) appear to have lower than aver- 
age masses (Figure [6J. These hosts also exhibit some of the 
largest ratios of M^ax/M^p, so the low masses may be at least 
partially due to their very young ages (4 of the 5 hosts are con- 
sistent with ages of only ~ 10 Myr; see Tableland Figure[U). 

The SSP stellar population ages span about 10 to 570 Myr, 
with a median value of (tssp) ~ 65 Myr. We find no clear 
correlation between the ages and redshifts of these hosts, al- 
though there appears to be a larger spread in ages at z < 
0.5 than at higher redshifts. The inferred ages are in good 
agreement with the analysis of IChristensen et aTl (l2004h for a 
smaller sample. 

Finally, we pay particular attention to the hosts of GRBs 
060505 and 060614 since these nearby bursts had long dura- 
tions but lacked an associ ated supernova, possibly pointing to 
a non - massive star origin (jDella Valle et al.|20 06; Fvnb o et ail 
2006 1 iGaFYam et alJ l200ff iGehrels et al J l200a lOfek et alJ 
2007[ ILevesque & Kewleyll2007t iThone et al]l2008l) . For the 
host of GRB 060505 we find a best-fit 3 SSP age of about 10 
Myr, with a stellar mass of about 2 x 10 8 M Q , while for the 
host of GRB 060614 we infer an age of 0.57 Gyr and a mass 
of 1.3 x 10 8 M Q (TableQ. The maximal masses are 5 x 10 9 
Mq and 1.6 x 10 9 M Q , respectively. The inferred masses ex- 
hibit better agreement with the distribution of long GRB hosts 
than with short GRB hosts. The very young stellar popu- 
lation age for the host of GRB 060505 is similar ly reminis- 
cent of the long GRB h ost population (see also lOfek et all 
I2007t lThone et alJl2008h . while the SSP age for the host of 
GRB 060614 is the largest of any of the long GRB hosts in 
this paper, and more closely in line with the short GRB host 
population. 

5. DISCUSSION 

Having established the basic stellar population properties 
of the short GRB hosts in our sample we now address the 
fundamental questions listed in |T] through a comparison to 
the general galaxy mass function and the hosts of long GRBs. 

5.1. Comparison to the Mass Distribution of Galaxies 

The question of whether the short GRB rate tracks stellar 
mass alone can be addressed through two statistical tests: (i) 
the fraction of short GRBs in early- and late-type galaxies; 
and (ii) the relation of the host stellar masses to the general 
galaxy mass function. 

3 This includes a rest-frame extinction of Ay.host =0.3 mag. A model with 
no extinction provides a poorer fit (%r = 2.4) and leads to an older age of 0.23 
Gyr and a larger mass of 2.5 X 10 9 Mq. 
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5.1.1. Host Demographics 

In the redshift range under consideration here (z ~ 0.2- 1) 
about an equal fraction of the cosmic stell ar mass budget re - 
sides in early- and late-type galaxies (e.g.. lllbert et al.l l2010). 
Therefore, if short GRBs track stellar mass alone we expect 
a roughly one-to-one ratio of galaxy types. This does not ap- 
pear to be the case. For example, within the sample of short 
GRBs with optical afterglows (19 events), only 2 are unam- 
biguous ly hosted by early -type galaxies (GRBs 050724 and 
1 00 1 1 7 ; iBerger et al Jl2005l Fong et al. in prep.), while 8 are 
unambiguously hosted by late-type galaxies; the probability 
of obtaining this ratio from an intrinsic one-to-one distribution 
is only 0.04. The identity of the remaining 9 hosts is unclear 
at the present due to their faintness or the lack of underlying 
galaxies at the burst positions. Still, unless nearly all of these 
bursts were hosted by early-type galaxies, the resulting ratio 
appears to be skewed in favor of late-type host galaxies with 
on-going star formation activity. We note that the same result 
holds true if we consider the bursts with only X-ray afterglow 
positions and identified hosts. 

Thus, the host galaxy demographics suggest that short 
GRBs do not track stellar mass alone, or phrased alternatively, 
they do not have a delay time distribution that is skewed to old 
ages of ~ few Gyr. It is possible, however, that this result is in- 
fluenced by secondary factors such as the typical circumburst 
density or intrinsic differences in the energy scale and after- 
glow brightness as a function of galaxy type (possibly remi- 
niscent of the differences in peak lumino sity for Type la su- 
pernovae in early- an d late-type galaxies; lHamuv et al.ll2000b 
iMannucci et alj[2006l) . If such differences lead to fainter af- 
terglows (or prompt emission) for short GRBs in early-type 
galaxies, this would suppress the early-type fraction. Al- 
though the modest size of the host sample, and the substantial 
fraction of short GRBs with only 7-ray positions (~ 1 /3 of all 
events), prevent definitive conclusions, it does not appear that 
the optical afterglo ws of short GR Bs in early- and late-type 
galaxies are distinct (Berger 2010a|). 

Thus, our preliminary conclusion from the host demograph- 
ics is that short GRBs do not appear to track the cosmic stellar 
mass fraction in early- and late-type galaxies. 

5.1.2. Comparison to the Galaxy Mass Function 

We next turn to a comparison of the inferred stellar masses 
with the galaxy mass function. The cumulative distribution 
of stellar masses for the short GRB hosts is shown in Fig- 
ure|5]with the range of possible masses bounded by the single 
age SSP and maximal values. We also present a breakdown 
of the sample into early- and late-type galaxies, each span- 
ning the same range. For the late-type hosts, the intermediate 
young+old values are also shown. To compare these distri- 
butions to the distribution of galaxy masses we also plot the 
cumulative stellar mass function weighted by mass, i.e., the 
fraction of stellar mass in galaxies above some mass, /(> M), 
given by the equation: 



/(>M) = 



f™M' x $(M')dM' 



/ °°M' x $(M')dM' 
where <£>(M) is the Schechter mass function: 



$(M) = $* 



M 

M* 



exp 



M 

M* 



(2) 



(3) 



2MASS/2dF catalogs for all galaxy types at z ~ (AT = 
10 UM M Q , a = -1.18); the nearly identical iPanter et ail 
(2004) mass function from SDSS for al l galaxy types a t 
z ~ (M* = 10 1119 M , a = -1.16); the iBell et 



/ t ypes a t 
1JO003) 



We use several determinat i ons of <&(M) in this compari- 
son, including the ICole et all d2001l) mass function from the 



mass function for late-type galaxies from 2MASS/SDSS con- 
verted to a Salpeter IMF for comparison with our inferred 
values (M* = 10 10 97 M , a = -1.27); and the lllbert et all 
(2010) mass functions from the COSMOS survey for quies- 
cent galaxies at z ~ 0.3 (M* = 10 11 13 M , a = -0.91) and 
intermediate-activity galaxies at z ~ 0.5 (M* = 10 10 93 M , 
a = -1.02), matched to the redshifts of the early- and late- 
type short GRB host galaxies in our sample. The result- 
ing distributions of /(> M) for the various mass functions 
are shown in Figure [5] The median of each distribution is 
clo se to the value o f M* , but also depen ds on the value of 
a (ICole etaimOQll : (AT) = 10' 88 M ; IPanter et al.ll200l 
(M*) = IP 1092 M ; IBell etalJl2003l (AT) = 10 10 61 M ; 
lllbert et all20Tol (AT*) = 10 10 76 M for intermediate-activity 
and (AT) = 10 11 02 M for quiescent). 

In $4] we demonstrated that the median single age SSP 
mass for the short GRB host sample is 1.3 x 10 ll) M , or 4 
ps 0.2(Af*), well below the expected value if the short GRB 
rate tracked stellar mass. The result is even more discrepant if 
we restrict the sample to only the late-type hosts, with a me- 
dian of about 0.06 (A/*). On the other hand, for the early-type 
hosts we find a median of about 1 .6(Af*), suggesting that they 
are drawn from the early-type galaxy mass function. 

Using the maximal masses for the full sample we find a 
median stellar mass of about 1.2(M*), while for the subset 
of late-type hosts the median is about 0.5 (Af*); the maximal 
masses of the early-type hosts lead to a large median mass 
of about 5 (A/*), suggesting that they are likely over-estimates 
of the true stellar mass. Finally, using the young+old stel- 
lar masses for the late-type hosts the median mass is about 
0.3(Af*). Thus, we find that the distribution of stellar masses 
for the late-type hosts is shifted to a somewhat lower value 
than the expected median of the galaxy mass function, while 
the early-type hosts have a similar median to the galaxy mass 
function. 

Beyond the comparison of median values, we further assess 
the agreement (or lack thereof) between the short GRB hosts 
and galaxy mass functions using the Kolmogorov-Smirnov 
(K-S) test. For the full sample there is negligible probabil- 
ity that the distribution of single age SSP masses is drawn 
from the galaxy mass function, with P«8x 10~ 5 . On the 
other hand, for the maximal mass distribution the probability 
is P ps 0.6 indicating that for these masses the short GRB sam- 
ple is fully consistent with the galaxy mass function. Using 
the intermediate case of SSP masses for the early-type hosts 
and the young+old masses for the late-type hosts, we find a 
probability of about 0.3, indicating that this combination is 
also fully consistent with the galaxy mass function. 

Separating the early-type hosts we find that thei r SSP 
masses are fully consistent with the lllbert et alJ {2010) mass 
function of quiescent galaxies (P ps 0.8); their large maxi- 
mal masses, on the other hand, are inconsistent with the mass 
function, with P ps 0.007. Finally, for the late-type hosts we 
fin d a clear inconsist ency of the single age SSP masses with 
the lIlbertetaT](l2010l) mass function of intermediate-activity 
galaxies, with P ps 4 x 10~ 7 . However, the maximal mass dis- 
tribution is fully consistent with the mass function (P ps 0.3), 

4 In the conversion to units of (M* ) we use the mass function appropriate 
to each galaxy type, or the combined sample. 
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while the young+old mass distribution is marginally consis- 
tent (P« 0.1). 

To summarize, the distribution of short GRB host masses 
is compatible with the overall mass distribution of galaxies 
only if their stellar masses are given by the SSP masses for 
the early-type hosts and the maximal or young+old masses 
for the late-type hosts. Since the opposite scenario (maximal 
masses for the early-type hosts and SSP masses for the late- 
type hosts) is unlikely, we conclude that the existing sample 
of short GRB hosts is consistent with the galaxy mass func- 
tion. Equivalently, this means that short GRBs may indeed 
track stellar mass alone. However, we caution that the host 
demographics seem to be at odds with the expected equal frac- 
tions of total stellar mass in early- and late-type galaxies, un- 
less nearly all of the unidentified hosts are early-type galaxies. 
This, along with the somewhat lower than expected masses of 
the late-type hosts, leaves open the possibility that at least a 
subset of short GRB progenitors track star formation activity 
rather than stellar mass. 

5.2. Comparison to Long GRBs 

Despite the possibility that some short GRB progenitors 
may track star formation activity, we find that the inferred 
stellar masses and population ages of short GRB hosts are 
generally distinct from those of long GRB hosts in both the 
single age SSP and maximal models. Most importantly, this is 
true for the subset of late-type hosts. In the framework of the 
single age SSP model we find a K-S probability of only 0.006 
that the long and short GRB hosts are drawn from the same 
mass distribution. The probability is higher for the subset of 
late-type short GRB hosts, PrjO.1 (Figure[3]i. However, since 
the SSP values represent the mass of only the young stellar 
populations, they are mostly reflective of the star formation 
activity rather than the total stellar mass. If we use instead 
the maximal masses, the K-S probability that the long GRB 
hosts and late-type hosts of short GRBs are drawn from the 
same sample is negligible, P w 4 x 10~ 5 (Figure [H, demon- 
strating that they are distinct galaxy populations. A similar 
conclusion is apparent from a comparison of the single age 
SSP population ages. The K-S probability that the long GRB 
hosts and late-type hosts of short GRBs are drawn from the 
same distribution is only P w 0.006. 

Thus, the long GRB hosts have significantly lower stel- 
lar masses than the subset of late-type short GRB hosts, and 
their young stellar population are significantly younger. In- 
de ed, a comparison o f the long GRB host maximal masses to 
fhe lllbertetal1(l2010h mass function of high-activity galaxies 
at z ~ 0.7 (appropriate for the long GRB sample considered 
here) indicates a K-S probability of only 0.002 that the long 
GRB hosts are drawn from the galaxy mass function. This 
is consistent with our understanding that their massive star 
progenitors select galaxies by star formation (and perhaps ad- 
ditional factors such as metallicity). 

The apparent distinction between long GRB hosts and the 
late-type hosts of short GRBs in terms of their stellar masses 
and young population ages strengthens our similar previous 
conclusion based on their star formation rate s, specific sta r 
formation rates, luminosit ies, and metallici ties ((Berger 2009), 
as well as physical sizes (Fon g et aT1l2010h . In essentially ev- 
ery property the late-type short GRB hosts point to a popu- 
lation of more quiescent, massive, and evolved galaxies than 
the hosts of long GRBs. We conclude that this rules out the 
idea that short GRB progenitors in late-ty pe hosts are massiv e 
stars identical to long GRB progenitors ( Virgili et al. 2009), 



even if the short GRBs in late-type galaxies indeed track star 
formation rather than stellar mass. 

5.3. The Delay Time Distribution 

A determination of the delay function from the derived stel- 
lar population ages is complicated by two primary factors. 
First, we have to assume that the short GRB progenitors in 
each host were formed within the inferred stellar population. 
This assumption is justified statistically both for an associa- 
tion of the progenitors with stellar mass and with star forma- 
tion activity, as long as we can appropriately normalize the 
rates of short GRBs. Second, while we can determine single 
age SSP ages from the broad-band photometry, these data are 
not sufficient to provide an age breakdown (by mass) for mul- 
tiple stellar components. Indeed, for our hybrid young+old 
model we had to fix the age of the old population (to the 
age of the universe, in this case). Still, we find that in the 
young+old model, the bulk of the mass (« 55-99%; Tabled 
is contained in the old stellar population, and so the progen- 
itors would have "old" ages (r > tssp) if they tracked stellar 
mass. 

As a result of these limitations we can only explore the 
implications of two main scenarios, namely that short GRBs 
track mass and/or star formation ac tivity . In the context of the 
former scenario we have shown in 35. ll that the short GRBs in 
early-type hosts trace stellar mass. Therefore, we can use their 
SSP ages to provide a rough estimate of the progenitor ages, 
which we find to be r w 0.8-4.4 Gyr, with a median of about 
3 Gyr. On the other hand, for the late-type hosts (for which 
we can extract no credible information on the mass-weighted 
stellar population age), we can infer a typical delay relative to 
the most recent star formation episode under the assumption 
that these progenitors track star formation activity. We find 
SSP ages of r w 0.03-0.5 Gyr, with a median of about 0.25 
Gyr, or young+old ages of about 0.01 -0.4 Gyr with a median 
of about 0. 1 Gyr. 

Thus, if short GRBs follow both stellar mass (in early-type 
galaxies) and star formation activity (in late-type galaxies), 
the typical delay times are about 3 and 0.2 Gyr, respectively. 

6. CONCLUSIONS 

Using multi-band photometry for 19 short GRB host galax- 
ies we extracted stellar masses and ages under the assumption 
of single age populations, combined young+old populations, 
and maximal mass populations. The resulting values allow 
us to investigate whether short GRBs track the cosmic stellar 
mass distribution, to estimate their typical delay times if mass 
and/or star formation play a role, and to investigate whether 
short GRBs in late-type galaxies share progenitors with long 
GRBs. The basic result is that the early-type host SSP masses 
track the mass distribution of early-type galaxies, while for 
the late-type hosts this is only the case if they have maxi- 
mal or possibly young+old masses. However, the host demo- 
graphics do not appear to follow the expected mass-weighted 
one-to-one ratio, unless nearly all of the current hosts of un- 
known type are early-type galaxies. We note that additional 
support for a short GRB rate tracking stellar mass is provided 
by high resoluti on Hubble Space Telescope (HST) imaging 
(Fong et al.| [2010h . which indicates that the locations of these 
bursts better track stellar mass (rest-frame optical light) than 
star formation (rest-frame UV light). 

We further find that the late-type hosts of short GRBs and 
the hosts of long GRBs are distinct in their stellar masses and 
young population ages, with the former having higher masses 
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by about an order of magnitude and older ages by about a 
factor of 4. This result is in good agreement with the previous 
conclusion that short GRB hosts have lower star formation 
rates and hig her metallicities compared to long GRB hosts 
(Berge3|2002). 

These observational results lead to several crucial conclu- 
sions regarding the nature of short GRB progenitors: 

• Short GRB progenitors are generally consistent with 
arising in old stellar populations that track the cosmic 
stellar mass distribution. 

• Despite this overall consistency, the dearth of early- 
type hosts, and the somewhat lower than expected 
masses of the late-type hosts indicate that short GRB 
progenitors in late-type galaxies may partially track star 
formation activity. 

• If short GRB progenitors track only stellar mass, then 
the typical stellar population ages of the early-type 
hosts indicate a typical progenitor age (or delay time) 
of ~ 3 Gyr (with a range of ~ 1 -5 Gyr). If the progen- 
itors in late-type hosts track star formation activity, the 
resulting typical delay is about 0.2 Gyr. 

• Even if short GRBs in late-type galaxies track star for- 
mation activity, their delay times are at least a few times 
longer than for long GRBs, indicating that they do not 
share the same progenitors. This is also supported by 
the substantially different mass distributions of long 
and short GRB hosts. 

The conclusion that short GRBs are consistent with track- 
ing stellar mass provides additional support for the compact 
object merger model, although it does not rule out other pro- 
genitors with substantial time delays (e.g., binary white dwarf 
mergers, or white dwarf and neutron star accretion-induced 
collapse). In the context of these various models, it is possi- 
ble that a short delay channel exists, but a substantial popula- 
tion with delays of onl y a few Myr or even tens of Myr (e.g., 
Belczynski et al. 2006) is not supported by the data. 

The clear distinction between short and long GRB hosts, 
and the inferred delay times, call into question the results of 
previous analyses claiming a substantial overlap between long 
and short GRB progenitors. For example, it demonstrates 
that bl unt statistical to ols such as the logA^-logS distribu- 
tion (IVirgili et al.ll2009h do not have sufficient time resolution 
to distinguish between truly young populations such as long 
GRB massive star progenitors, and intermediate-age progeni- 
tors as inferred here if short GRBs in late-type galaxies track 
star formation activity. Similarly, these results call into ques- 
tion a recent proposal that some short GRBs are due to widely 
off-ax is "cocoon" emission from long GRBs dLazzati et al.l 
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Our results also demonstrate that caution should be taken 
with the proposed re-classification of short and long GRBs 
into Type I a nd II events, marking old and young progenitors, 
respectively (Zhang et al. 20071 [2009b . Such a new bimodal 
classification may lead to the erroneous conclusion that short 
GRBs in late-type galaxies (even if they track the on-going 
star formati on activity) share the same progenitors as long 
GRBs (e.g- IVirgili et~all 2009) since both would be classified 
as Type II. At the very least, such a new classification scheme 
may require a further breakdown of the Type II events into 
those that result from massive stars versus those that simply 
track star formation activity with a modest delay, e.g., Type 
Ha and lib. Clearly, this is beyond the scope of the current 
short GRB sample. We note, by analogy, that for Type la 
supernovae a connection with both ste llar mass and star for- 
mation activity has b een demonstra ted dMannucci et"aTll2006t 
ISullivan et alj|2006t lAubourg et al.ll2008tlBrandt et alj|2010h . 
but the prompt channel is not related to the massive star pro- 
genitors of Type II and Ib/c supernovae. 

We conclude by noting that two primary refinements of 
our analysis technique can be made with future observations. 
First, high signal-to-noise ratio spectra can be used to ex- 
tract more detailed star formation histories than the available 
broad-band photometry. Such reconstructed histories will al- 
low us to more robustly estimate the delay ti me distribution 
in the context of a mass-weighted rate (c.f., Aubourg et al. 
120081 for Type la supernovae). Second, high angular reso- 
lution multi-band imaging with HST will allow us to deter- 
mine spatially-resolved ages at the locations of short GRBs 
and hence to further refine the connection with young or old 
populations. Both of these techniques are within our grasp in 
the near future. 
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TABLE 1 

Short GRB Host Galaxy Photometry 
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060801 14 h 12 ra 01.35 s +16°58'53.7" 1.130 N 



061006 07 h 24 ra 07.66 s -79°11'55.1" 0.438 Y 



061210 09 h 38 m 05.27 s +15°37'17.3" 0.410 N 



061217 10 h 41 m 39.32 s 



-21°07'22.1" 



070429b 21 h 52 ra 03.84 s -38°49'42.4" 



0.827 N 



0.902 N 



070714b 03 B 51 m 22.30 s +28°17'50.8" 0.923 Y 



070724 01 h 51 m 13.96 s 



-18°35'40.1" 



0.457 Y 



SDSS 

SDSS 

SDSS 

SDSS 

SDSS 

2MASS 

2MASS 

2MASS 

Hiorth et al. (2005) 
Covino et al. (2006) 
Covino et al. (2006) 
Covino et al. (2006) 
This Work 
This Work 

Gorosabel et al. (2006) 
Gorosabel et al. (2006) 
Gorosabel et al. (2006) 
Gorosabel et al. (2006) 
Berger et al. (2005) 
Gorosabel et al. (2006) 
Gorosabel et al. (2006) 
Gorosabel et al. (2006) 
This Work 
This Work 
This Work 
This Work 
This Work 
his Work 

Soderbe rg et al. (2006) 

Soderberg et al. (2006) 

Soderbe rg et al. (2006) 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

D'Avanzo et al. (2009) 
D'Avanzo et al. (2009) 
This Work 

D'Avanzo et al. (2009) 
D'Avanzo et al. (2009) 
This Work 

D'Avanzo et al. (2009) 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

This Work 

Graha m et al. (2009) 

Graham et al. (2009) 

Graham et al. (2009) 

Graham et al. (2009) 

Graham et al. (2009) 

Graham et al. (2009) 

Graham et al. (2009) 

This Work 

This Work 

This Work 
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TABLE 1 — Continued 



GRB 


RA 


Dec 


2 


OA? 


Type a 


Filter 


m AB b 


Refs. 




(J2000) 


(J2000) 










(mag) 
















z 


20.26 ±0.04 


This Work 














J 


20.01 ±0.02 


This Work 














H 


19.78 ±0.02 


This Work 














K 


19.71 ±0.04 


This Work 


070729 


03 h 45 m 16.04 s 


-39° 19' 19.9" 




N 


? 


g 


24.38 ±0.38 


This Work 














r 


23.31 ±0.03 


This Work 














i 


21. 80 ±0.06 


This Work 














z 


21.82±0.17 


This Work 














J 


20.85 ±0.10 


This Work 














K 


20.10±0.10 


This Work 


070809 r 


13 h 35 m 04.4P 


-22°08'28.9" 


0.473 


Y 


E 


g 


21.80±0.05 


This Work 














r 


19.89 ±0.02 


This Work 














i 


19.27 ±0.05 


This Work 














K 


17.96 ±0.04 


This Work 


071227 


03 h 52 m 3 1.26 s 


-55°59'03.5" 


0.381 


Y 


L 


g 


22.82±0.13 


This Work 














r 


20.60 ±0.05 


This Work 














i 


20.47 ± 0.04 


This Work 














z 


19.77 ±0.03 


This Work 














J 


19.16±0.06 


This Work 














K 


18.15 ±0.06 


This Work 


080123 


22 h 35 ra 46.10 s 


-64°54'03.2" 


0.495 


N 


L 


g 


22.06 ±0.06 


This Work 














r 


20.89 ±0.05 


This Work 














i 


20.49 ±0.07 


This Work 














z 


20. 12 ±0.20 


This Work 














J 


20.30 ±0.05 


This Work 














K 


19.58 ±0.06 


This Work 


090510 


22 h 14 m 12.56 s 


-26° 34' 59.0" 


0.903 


Y 


L 


g 


23.78 ±0.08 


This Work 














i 


22.41 ±0.14 


This Work 














z 


22.66±0.17 


This Work 














J 


21.79±0.15 


This Work 














K 


21.98 ±0.35 


This Work 


090515 c 


10 h 56 ra 36.11 s 


+14°26'30.3" 


0.403 


Y 


E 


g 


21.89 ±0.02 


This Work 














r 


Z\).Zl ± U.Uj 


Ihis WOIK 














i 


19.45 ±0.05 


This Work 














J 


18.67 ±0.05 


This Work 














K 


18.23 ±0.05 


This Work 


100117 


00 h 45 ra 04.7 s 


-01°35'42.0" 


0.920 


Y 


E 


g 


26.01 ±0.30 


Fong et al. in prep. 














r 


23.97 ±0.04 


Fong et al. in prep. 














J 


21.72 ±0.26 


Fong et al. in prep. 














H 


21.63 ±0.21 


Fong et al. in prep. 














K 


21.24±0.25 


Fong et al. in prep. 



NOTE. — Properties and photometric measurements for short GRB host galaxies in this paper. 
Type indicates whether a galaxy is early-type (E) or late-type (L). 

Corrected for Galactic extinction iSchlegel et al. 1998). 

These bursts do not have coincident host galaxies. Magnitudes are provided for the galaxy with the lowest probability of chance coincidence (Berger 2010a). 
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TABLE 2 

Stellar Ages and Masses of Short GRB Host Galaxies 



GRB 


z 


TvDe 




SSP" 






Young+Old 


h 




loH(M„ rl .) f ' 


T 

(G\r) 


log(M) 
!Mn) 




Ty 

(Gvr) 


log(M Y ) 
(M ) 


(Gyr) 


lo2(Mn) 
(Mn) 


X 2 


050509b 


0.225 


E 


3.18 


11.6 


0.7 












11.9 


050709 


0.161 


L 


0.26 


8.8 


0.6 


0.21 


8.7 


11.6 


9.2 


0.6 


9.7 


050724 


0.257 


E 


0.94 


10.8 


0.3 












11.7 


051210 


1.3±0.3 e 


7 


0.03 


8.8 


2.5 














051221 


0.547 


L 


0.17 


9.4 


2.7 


0.14 


9.3 


8.3 


<9.4 


2.4 


10.1 


060801 


1.130 


L 


0.03 


9.1 


1.2 


0.03 


9. 1 


5.4 


< 10.8 


1.9 


10.9 


061006 


0.438 


L 


0.24 


9.0 


1.3 


0.09 


8.5 


O 1 

9. 1 


9.7 


0.8 


9.8 


061210 


0.410 


L 


0.38 


9.6 


4.2 


0.38 


9.5 


9.1 


< 9.9 


5.2 


10.5 


061217 


0.827 


L 


0.03 


9.1 


1.7 


0.03 


9.0 


6.6 


< 9.3 


2.1 


10.1 


U /U4ZyD 


yj.yvZ 


I 

L 


U.40 


1U.4 


J .ft 


U.J 3 


10.2 


6.3 


> 11 A 

< ll.U 


Z.O 


11. J 


070714b 


0.923 


L 


0.22 


9.4 


2.9 


0.10 


9.1 


6.2 


< 10.2 


3.3 


10.2 


070724 


0.457 


L 


0.30 


10.1 


2.1 


0.01 


8.3 


8.9 


10.9 


0.3 


11.0 


070729 


0.8±0.1" 


? 


0.98 


10.6 


2.7 














070809 


0.473 


E 


3.07 


11.4 


0.9 












11.7 


071227/ 


0.381 


L 


0.49 


10.4 


2.2 


0.36 


10.3 


9.3 


< 11.0 


2.7 


11.4 


080123 


0.495 


L 


0.31 


10.1 


1.7 


0.03 


9.7 


8.7 


< 10.7 


1.8 


11.1 


090510 


0.903 


L 


0.14 


9.7 


3.3 


0.02 


9.3 


6.3 


< 10.4 


3.7 


10.6 


090515 


0.403 


E 


4.35 


11.2 


0.1 












11.5 


100117 


0.920 


E 


0.79 


10.3 


1.8 












10.9 



NOTE. — Stellar population parameters derived from the broad-band SEDs in Table[T]using the Maraston 12005) stellar population synthesis models (^3). 
" Fit with a single stellar population age. 

* Fit with a two-component young+old SSP model, with the old population fixed at an age of to equivalent to the age of the universe at each host redshift, and 

a freely- varying young population with an age of Ty. 

c to is fixed at the age of the universe at the appropriate redshift. 

d The maximal mass is derived using the observed A"-band fluxes assuming a stellar popul ation age equal to the age of the universe at the host redshift. 

' A spectroscopic redshift is not known; these values are derived photometrically from the Maraston (2005) models. 

' This model includes host extinction with Ay . host = 1-5 mag; a model with no host extinction does not provide a satisfactory fit. 



12 



Leibler & Berger 



TABLE 3 

Stellar Ages and Masses of Long GRB Host Galaxies 



SSP 



GRB 


z 


Age 


Mass 


2 

XT 


log(M ma x)* 






(Gyr) 


(M ) 




(Mq) 


970228 


0.695 


0.14 


8.6 


1.1 


8.7 


970508 


0.835 


0.07 


8.5 


0.4 


8.8 


98061 3 


1.097 


0.05 


9 j 


3.6 


9.4 






0.02 c 


9.0 


1 .0 




980703 


0.966 


0.05 


9.3 


1.0 


10.8 


990123 


1.600 


0.05 


9.4 


2.0 


9.7 


990712 


0.433 


0.11 


9. 1 


1.1 


9.3 


991208 


0.706 


0.10 


8.8 


0.6 


9.0 


000? 1 


0.846 


0.12 


9.3 


0.7 


9.5 


00041 8 


1.118 


0.05 


9.3 


2.0 


9.6 


0009 1 1 


1.058 


0.06 


8.9 


0.2 


9.6 


010222 


1.480 


0.06 


8.6 


0.8 


8.9 


010921 


0.451 


0.22 


9.6 


0.4 


9.6 


011121 


0.362 


0.24 


9.5 


2.0 


10.3 


020813 


1.255 


0.09 


9.4 


1.3 


9.6 


020819b 


0.410 


0.35 


10.3 


1.0 


10.5 


020903 


0.251 


0.01 


7.8 


0.9 


9.1 


030328 


1.520 


0.06 


9.3 


0.7 


9.6 


030329 


0.168 


0.01 


6.8 


1.1 


8.0 


030528 


0.782 


0.01 


9.0 


0.4 


10.0 






0.04 


8.6 


0.7 




050826 


0.296 


0.50 


9.9 


2.0 


11.0 


060218 


0.034 


0.01 


6.7 


0.1 


9.1 






0.06 


7.5 


0.3 




060505 


0.089 


0.23 


9.4 


2.4 


9.7 






0.01 d 


8.3 


0.9 




060614 


0.125 


0.57 


8.1 


0.3 


9.2 



NOTE. — Stellar population parameters derived from the broad-band SEDs using the Maraston 12005) stellar population synthesis models (^3j. 
" Fit with a single stellar population age. 

* The maximal mass is derived using the observed A'-band flux in conjunction with the maximal age allowed by the burst redshift. 
c This model includes host extinction with Ay host = 0.5 mag. 
d This model includes host extinction with Ay.host =0.3 mag. 
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FIG. 1 . — Optical and near-IR spectral energy distributions of short GRB host galaxies (black circles). Each SED is fitted with a Maraston (2005) single stellar 
population model (red line) through a maximum likelihood fit of the synthesized photometry (red squares). The resulting best-fit mass and age are listed in 
Tabled 
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FIG. 2. — Optical and near-IR spectral energy distribution of GRB 050709 with the three models used in this paper to extract the stellar mass and population 
age. Symbols are as in Figure[T] Left: Single age SSP model. Center: Young+Old SSP model (magenta=old; cyan=young) with the old population age fixed at 
the age of the universe at the redshift of the burst (z = 0.161 in this case). Right: Maximal mass model with the population age fixed at the age of the universe 
and using only the A'-band photometry. The Young+Old model leads to total masses intermediate between the single age SSP and the maximal models, and has 
younger ages for the young population than the single age SSP model. The resulting best-fit masses and ages for the three models are listed in Table[2] 
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log [M / M ] 

a L SUIT 

FIG. 3. — Histograms of inferred stellar masses from the single stellar population fits shown in Figure[T]for the hosts of short (black) and long (gray) GRBs. 
The inset shows the cumulative distributions, including for the subset of late-type short GRB hosts (blue). The median values for the three samples are given 
in parentheses, and the Kolmogorov-Smirnov probabilities that the distributions of short and long GRB hosts, as well as star forming short GRB and long GRB 
hosts are drawn from the same distribution are provided in the inset. 
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FIG. 4. — Same as Figure[3]but for the maximal masses. 
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FIG. 5. — Cumulative distributions for the full sample of single age simple stellar population (SSP) masses, maximal masses, and combined young+old and 
SSP masses for the late- and early-type hosts, respectively (black; bottom panel). The upper panel shows a breakdown by galaxy type (late-type: blue; early- 
type: red). The shaded regions represent the range of possible stellar masses since the SSP masses, which are effectively light-weighted values, are most likely 
an under-estimate, while the maximal masses make the extreme assumption that all hosts are dominated by populations with the age of the universe. For the 
late-type hosts we also plot the total masses from a young +old SSP fit (TableEJ, which are more closely representative of the total mass. Also shown are the 
fractions of total stellar mass in galaxies with mass , > M, cal culated from se veral published gal axy stellar mass functions at z ~ 0-2 (cyan and magenta lines; 
ICole et al. 2001; Bell et al. 2003; Panter et al. 2004; Ilb ert et al.1201 J) ; for the lllbert et alj f2010t) mass function we use the z ~ 0.5 bin, appropriate for the short 
GRB sample (Table^, and separately plot the mass function for quiescent galaxies and for intermediate-activity galaxies, which resemble the intermediate star 
formation activity in short GRB hosts i Berger 2009). The comparison indicates that short GRBs trace galaxy mass only if the bulk of the late-type hosts have 
close to maximal masses. The subset of early-type hosts appears to faithfully trace the mass function of galaxies for the SSP-derived masses. 
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FIG. 6. — Single stellar population masses plotted as function of redshift for the hosts of short (black) and long (gray) GRBs. Circles and stars designate early- 
and late-type short GRB hosts, respectively. Filled symbols designate hosts identified through optical afterglow positions, thick open symbols designate hosts 
identified through coincidence with Swift/XKT error circles, and thin open symbols designate galaxies identified as potential hosts for short GRBs with optical 
afterglows based on chance coincidence probabilities (Berger 2010a). Open gray symbols designate the ambiguous long GRBs 060505 and 060614. In all cases 
the lines indicate the maximal masses. 
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log [age / Gyr] 

FIG. 7. — Histograms of inferred stellar population ages from the single stellar population fits shown in Figure[T]for the hosts of short (black) and long (gray) 
GRBs. The inset shows the cumulative distributions, including for the subset of late-type short GRB hosts (blue). The median values for the three samples are 
given in parentheses, and the Kolmogorov-Smirnov probabilities that the distributions of short and long GRB hosts, as well as star forming short GRB and long 
GRB hosts are drawn from the same distribution are provided in the inset. 



20 



Leibler & Berger 



10 l 



CD 

g)10" 1 

< 



10" 



^ o 
o o 



* V 



• Short GRBs (early) 
^ Short GRBs (late) 
■ Long GRBs 



it 

★ 
★ 



it it 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

Redshift 

FIG. 8. — Single stellar population ages plotted as function of redshift for the hosts of short (black) and long (gray) GRBs. Symbols are as in Figure[6] 



